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CN 1 Abstract 

i-C ■ We investigate the effect of early chemical freeze-out on radial and el- 

ry ^ liptic flow by using a fully three dimensional hydrodynamic model. We 

find that the time evolution of temperature and the thermal freeze-out 
temperature dependence of average radial flow are different from the 
results by using a conventional hydrodynamic model in which chemical 
equilibrium is always assumed. We also analyse the pt spectrum and 
£T) ' V2{pt) at the RHIC energy and consistently reproduce experimental data 

© : by choosing the thermal freeze " out temperature Tth = 140 MeV ' 
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^ ■ 1 Introduction 



In the conventional hydrodynamic models for analyses of relativistic heavy-ion 
collisions, both chemical and thermal equilibrium are assumed and the chem- 
ical freeze-out temperature T c h is the same as the thermal freeze-out temper- 
ature T t h which is so chosen as to reproduce the slope of p t (or m t ) spectra. 
For analyses of the pt slope at the RHIC energy by using the conventional hy- 
drodynamics, see, e.g., Refs. [1-5]. On the other hand, various analyses based 
on thermal models indicate the following picture of time evolution of hot and 
dense hadronic matter |J: the system first undergoes the chemical freeze-out 
at T c h where the observed particle ratios are fixed, continues to cool down and 
finally goes through the thermal freeze-out where the momentum distribution 
is fixed.Q We take account of this picture in a hydrodynamic model and discuss 
how the early chemical freeze-out affects the space-time evolution of fluids. We 
also analyse the particle spectrum at the RHIC energy. In the present work, 
we assume [ib = for simplicity. Although we neglect finite baryon effects, 
our hydrodynamic simulations are performed in a genuine three dimensional 
space. As a future work, we will include finite /xb in our model and analyse 
particle ratios and spectra at various collision energies [0, §]. 

^^The observed momentum distribution can be affected by feeding from resonance decay 
after thermal freeze-out. 
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2 Equations of State 

We construct three models equation of state (EOS) to discuss the space-time 
evolution of fluids. These models describe the first order phase transition 
between the QGP phase and the hadron phase at T c = 170 MeV. We suppose 
the QGP phase is composed of massless u, d, s quarks and gluons and that it 
is common to three models. For the hadron phase, we choose three different 
models EOS as follows. 

The first model is a conventional resonance gas model in which complete 
chemical equilibrium is always assumed (model CE). We include hadrons up 
to the mass of A(1232) for simplicity. 

The second model is the simplest one which describes the picture of chemi- 
cal freeze-out (model CFO). Below T c h, we assume the numbers of all hadrons 
Ni included in the EOS are fixed and that the particle number densities obey 
<9 M nf = 0. We introduce a chemical potential /ii(T) associated with each 
species so that Ni becomes a conserved quantity. From the conservation of en- 
tropy <9 M s M = 0, the ratio of the particle number density to the entropy density 
obeys 



rii(T < T ch ) rn(T ch ) 



s(T < T c 



chJ 



s(T t 



(1) 



ch) 



for all hadrons along the adiabatic path. Here we assume T c h = 170 MeV, 
which is consistent with a recent analysis based on a thermal model || . From 
Eq. ([]]), we obtain a chemical potential as a function of temperature for each 
hadron. 
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Fig. 1: (Left) Pressure as functions of energy density for three models (Right) Temperature 
as functions of energy density. 

The third model represents a more realistic EOS than the second one. The 
observed particle numbers are always composed of the direct particles and the 
contribution from resonance decays, e.g., N n — N^ + di^N. So some 
processes with large cross sections can be equilibrated even below T c h as long 
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as the equality fn{T < T c h)/s(T < T c h) = nj(T c h)/s(T c h) is kept instead of 
Eq. (H). We assume 7r, i\T, 77, N, A and £ are "stable" particles and that all 
chemical potentials can be represented by chemical potentials associated with 
these stable particles, e.g., \i p = 2fi n , fix* — ^ + Hk, A*a = fJnr + /^tv an d so on 
[ prOfl . Thus the third model describes the partial chemical equilibrium (model 
PCE) even below T ch . 

One can easily evaluate EOS's for these models by using chemical potentials 
as obtained above. We represent pressure and temperature as functions of 
energy density for three models in Fig. 1. We find pressure as a function 
of energy density is similar to each other. On the other hand, temperature 
as a function of energy density in the models CFO or PCE is deviated from 
the model CE. See also Fig. 1 (right). Since the resonance population of the 
models CFO or PCE is larger than that of the model CE due to the chemical 
freeze-out, the energy density at a fixed temperature in the hadron phase is 
also large in those models. Conversely, temperature in the models CFO or 
PCE at a fixed energy density is smaller than in the model CE. 



3 Radial and Elliptic Flow 

By using the EOS's in the previous section, we perform the hydrodynamic 
simulation with appropriate initial conditions and compare the hydrodynamic 
evolution with each other. We choose initial conditions so as to roughly re- 
produce pseudorapidity distribution in Au+Au 130A GeV central (0-6 %) and 
semi-central (35-45 %) collisions obtained by the PHOBOS Collaboration [11 . 
At the initial time tq = 0.6 fm/c, the energy density (temperature) for cen- 
tral collisions at the centre of the fluid, where it has a maximum value, is 
30.8 GeV/fm 3 (349.5 MeV). For further details of our hydrodynamic model, 
especially initial longitudinal and transverse profiles, see Ref. PJ. 
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Fig. 2: (Left) Time evolution of temperature for three models. (Right) Thermal freeze-out 
temperature dependence of average radial flow for three models. 
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Figure 2 (left) shows the time evolution of temperature at the centre of 
the fluid for three models. From EOS's represented in Fig. 1 (left), we expect 
the time evolution of energy density is similar to each other. However, when 
we transpose the energy density to the temperature by using EOS's in Fig. 1 
(right), the difference between three models appears in the time evolution of 
temperature. We find temperature in the models CFO or PCE cools down 
more rapidly than in the conventional model CE. We next show in Fig. 2 
(right) the thermal freeze-out temperature dependence of average radial flow 
at vanishing space-time rapidity r] s = and compare our results with the 
ones from thermal models. Radial flows in the models CFO and PCE are 
considerably reduced. This is understood qualitatively by the same reason as 
the above discussion on the time evolution of temperature. These results seem 
to be consistent with the one obtained by Burward-Hoy [12], while these are 
completely inconsistent with the model II in Ref. flTBf . Since the results from 
thermal models are different from each other and we do not take into account 
the baryonic chemical potential, more detailed analyses based on both thermal 
and hydrodynamical models are needed to understand the thermal freeze-out. 
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Fig. 3: Thermal freeze-out temperature dependence of transverse momentum spectrum for 
negative charged hadrons (left) and V2{pt) for charged hadrons (right). Solid, dashed and 
dotted lines correspond to, respectively, Tth = 140, 120 and 100 MeV 

Finally we show in Fig. 3 the p t spectrum and t^G *) only for the model 
PCE. Left figure represents the transverse momentum distribution of negative 
charged hadrons for various thermal freeze-out temperatures. Experimental 
data is observed by the STAR Collaboration (jTJJ] . Although the p t slope usually 
depends on T t h in the conventional EOS where the chemical equilibrium is 
assumed, our result for the model PCE is found to be almost independent of 
Tth due to the reduction of radial flow. Figure 3 (right) shows the transverse 
momentum dependence of v-i for charged hadrons in minimum bias events. 
Our result with T t h = 140 MeV gives the best fit to the data obtained by the 
STAR Collaboration |l5j below 1 GeV/c. The discrepancy of T t h between our 



result and the one in Ref. [E| may be due to the fact that the resonance decays 
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and the chemical potential of pions are not taken into account in Ref. [|1^]. We 
consistently reproduce experimental data of both p t spectrum and v 2 (pt) below 
1 GeV/c by choosing T t h = 140 MeV, so contribution from hard processes can 
become dominated above 1 GeV/ c. Note that the concavity of V2{pt) in small p% 
region (p t < 0.3 GeV/c) is due to the resonance decay after thermal freeze-out 
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4 Summary and Discussion 

We have investigated the effect of early chemical freeze-out on radial and el- 
liptic flow at the RHIC energy We found the system cools down more rapidly 
than the conventional model due to the large population of resonances. We 
also found that the resultant radial flow is considerably reduced when we com- 
pare with each other at a fixed Tth. Although the space-time evolution of 
temperature in the model PCE is considerably different from the conventional 
model, we consistently reproduced p t spectrum and v 2 (pt) below 1 GeV/c by 
choosing the thermal freeze-out temperature T t h = 140 MeV. This indicates 
the contribution from hard processes can become dominated above 1 GeV/c. 
It should be noted that the average radial flow results in < v r >|^ s= o~ 0.35c. 

Since the two pion correlation function is believed to be sensitive to the size 
and the lifetime of the fluid, it is interesting to see the effect of early chemical 
freeze-out on the HBT radii. Early chemical freeze-out in hydrodynamics will 
be one of the possible explanations to interpret the "HBT puzzle" [[|] at the 
RHIC energy |17|, |j~8f . In addition to the HBT radii, detailed analyses of 



particle ratio and spectra with finite /xb are also our future works. 
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